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Abstract 



We analytically calculate higher order corrections to coefficient functions of the operator 
product expansion (OPE) for the Euclidean correlator of two energy-momentum tensors in 
massless QCD. These are the three- loop contribution to the coefficient Co in front of the 
unity operator Oq = 1 and the one and twodoop contributions to the coefficient C\ in front 
of the gluon "condensate" operator 0\ = — jG^G^. For the correlator of two operators 0\ 
we present the coefficient C\ at twodoop level (the coefficient function Cq is known at four 
^ ! loops from []]]). 



> '. 1 Motivation 

*n | The energy-momentum tensor correlator 

T^'P°(q) = i Jd 4 xe iqx (0|f^ CT (2;)|0), f^ pa {x) = T[T^ \x)T pa '(0)] (1) 

plays an important role in many physical problems. A lot of these lie in the field of Quark 
Gluon Plasma (QGP) physics. Here the correlator eq. ([T]) is the central object for describing 
transport properties, like the shear viscosity of the plasma (see e.g. 0, 01) and spectral 
functions for some tensor channels in the QGP fl. Another application is a sum rule approach 
to tensor glueballs. These special hadrons without valence quarks are determined by their 
gluonic degrees of freedom. QCD allows for such particles but a conclusive discovery has not 
yet been made. 

In a sum rule approach 0] one usually starts with the vacuum correlator of an interpolating 
local operator which has the same quantum numbers as the hadrons we want to investigate. 
If we are interested in glueballs we take local operators consisting of gluon fields. For the 
cases J PC = ++ ,0 ^ and 2 ++ the following operators are usually considered: 

Oxix) = —G^G^x) (scalar) (2) 

O x (x) = G^G^x) (pseudoscalar) (3) 

O t {x) = T^ u {x) (tensor) (4) 
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where is the gluon field strength tensor and 

G^i/ = £^vpuG p (5) 

the dual gluon field strength tensor. For more details see [f|. The vacuum expectation value 
(VEV) of the correlator of such a local operator 0(x) 

U(Q 2 ) = ijd 4 x e iqx (0|T[ O(x)O(0)} |0) (Q 2 = -q 2 ) (6) 

can of course be calculated in perturbation theory for large Euclidean momenta, but this is 
not enough. Starting from the perturbative region of momentum space we can probe into 
the non-perturbative region by means of an OPE. The idea originally formulated in 0] is to 
expand the non-local operator product i j d 4 x e tqx T [O(x)O(0)] in a series of local operators 
with Wilson coefficients depending on the large Euclidean momentum q. In sum rules we 
usually have dispersion relations connecting the VEV of such a Euclidean operator product 
to some spectral density in the physical region of momentum space. As we are ultimately 
interested in the VEV of this operator product we only have to consider gauge invariant 
scalar operators in the expansion. 

Effectively this expansion separates the high energy physics, which is contained in the Wilson 
coefficients, from the low energy physics which is taken into account by the VEVs of the local 
operators, the so-called condensates jjj. These cannot be calculated in perturbation theory, 
but need to be derived from low energy theorems or be calculated on the lattice. Such 
an OPE has already been done for the cases eq. ([2} and eq. ([2D (see @, @]) with one-loop 
accuracy. 

In this work we present the results for the Wilson coefficients in front of the operators Oo 
and [0\] for the correlator eq. ([T]) in massless QCD: 

TH"->f»( q ) = C^ pa {q)t + Cr ,pa (q)[Oi] + .... (7) 

q 2 — > — oo 



The brackets in [0\] indicate that we take a renormalized form of the operator 0\: 

[ 0l ] = Z G Of = -^G B ^G% (8) 

where the index B marks bare quantities. We start our calculation with bare quantities which 
are expressed through renormalized ones in the end: 

T^(q) = Cf^ pa (q)Of = Cr pa (q)[Oi}. (9) 

i i 

All physical matrix elements of [0±] are finite and so is the renormalized coefficien10 

a = ±c?. (io) 

The renormalization constant 

? s = l + .*ln^ = (l-^)" 1 (11) 



1 This statement as well as eq. (J9j) are only true modulo so-called contact terms; see a detailed discussion 
in the next section. 
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has been derived in a simple way in [nil ] (see also an earlier work 11]). Here Z aa is the 
renormalization constant for a s and we defin^l 



2 d 
p — T lna s 
d/x^ 



i>0 



i+1 



(12) 



In the massive case the two and four dimensional operators 0/ = and 2 = m ftp ftp f 
would have to be included as well for every massive quark flavour f. The VEVs of all other 
linearly independent scalar operators of dimension four vanish either by some equation of 
motion or as they are not gauge invariant. The contributions of higher dimensional operators 
are suppressed by higher powers of in the coefficients. 

Apart from the leading coefficient Cq in front of the local operator 0$ = 1 the coefficient 
C\ in front of [0\] is of special interest for many applications. One example is if we have a 
spectral density defined by our correlator and we want to calculate the shift in this spectral 
density from zero to finite temperature: 

Ap(u,T) = p(u,T) - p(u,0). (13) 

The spectral density at T = is calculated from the VEV of the correlator whereas for the 
spectral density at finite T we take the thermal average of the operator product. For the 
unity operator Oq = 1 the VEV and the thermal average are both 1 due to the normalization 
conditions. Hence the leading term from the OPE, i.e. the one proportional to Oq vanishes in 
eq. (|13|) which makes the Wilson coefficients in front of 0\ and 2 the leading high frequency 
contributions to eq. (|13|) . For more details see e.g. [121 ] . 



2 The energy- momentum tensor in QCD 



The energy-momentum tensor which can be derived from the Lagrangian of a field theory is 
an interesting object by itself. To be identified with the physical object known from classical 
physics and general relativity it has to be symmetric as well as conserved. A very general 
method to derive such an energy- momentum tensor can be found e.g. in This has 



firstly been done for QCD in [16(] and the result derived from the renormalized Lagrangian 



C 



1 



1 



z 3 g^ct - — (d^y + z 3 d p cdp c + g a z x d p c {A" x c ) 



+ % -Z 2 iP i(> + g.ZwtfjTip 



2A 

hip 



(14) 



2 Often in the literature Z as is used instead of Zg and OLsG^G^u instead of 0\. This is justified because 
up to first order in a s the renormalization constants Zg and Z as are the same. Only in higher orders Zg and 
Z aB differ and therefore Zg has to be used in such cases. 
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IS 



T pu = - Z 3 G W GP + ^d p AP)A u + ±(d u d p AP)A p 

+ Z 3 (d p cd u c + d v cd^c) + gZ\ (d p c(A u x c) + d u c(A p x c)) 

+ (K^v + o^7m) ^ + (^m T 7^ + A v Tj p ) if; ^ 



9tiu 



|-iz 3 G pa G pa + i {d a d p A") A* + ^ 



Here 



+ Z 3 ^c^c + 5 Zi (3 p c(A p x c)) + % -Z 2 ^t^ + gZ x ^4^^\ 



G pu = d p A u - d v A p + (Ap xA u ), (16) 

Z3 



Z 3 ,Z 3 and Z2 stand for the field renormalization constants for the gluon, ghost and quark 
fields respectively and Z\ and Z\^ for the vertex renormalization constants. The abbreviation 
(A p x A v ) a = f abc A b ^Al, where f abc is the structure constant of the SU(iV c ) gauge group, is 
used and all colour indices are suppressed for convenience. 

This energy-momentum tensor consists of gauge invariant as well as gauge and ghost terms. 
If we were to consider general matrix elements of operator products we would have to include 
all these terms. It has been pointed out in 16[ however that for Green's functions with only 
gauge invariant operators it would be enough to take the gauge invariant part of the energy 
momentum tensor: 

T p i* I g i nv ZzGppG? + ^Z 2 ip (<^!> + 0^7/^) tp + %Z\^) {A p T-/„ + A v Tj^) -0 

(17) 

g^v 1 



X -Z z G pa G pa + l -Z 2 ^$^ + gZ^ATxjA . 



This has been checked in our calculation of Co which we have done once with the full energy- 
momentum tensor eq. (fT5j) and once with the gauge invariant part eq. (fT?]) up to three-loop 
accuracy. As expected both calculations yield the same result. 

The insertion of a local operator into a Green's function corresponds to an additional vertex 
in every possible Feynman diagram. For the energy-momentum tensor eq. (|15|) we get the 
vertices shown in Figure [Q 




Figure 1: Energy- momentum tensor vertices and there dependence on g s 
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In [16| it has been proven that the energy-momentum tensor of QCD is a finite opera- 
tor which means the Z-factors appearing in (|15p make any Green function of (renormal- 
ized) QCD elementary fields with one insertion of the operator T^ IU finite. We have 
used this theorem as a check for our setup and have calculated one - and two-loop cor- 
rections to the matrix elements (gluon(p,//i)|T^|gluon(p,/i2)), (ghost(p)|T Atl/ |ghost(p)} and 
(quark(p), gluon(0,^i)|T^|quark(p)) which turned out to be finite as expected. 

Another important consequence of the finiteness property is the absence of the anomalous 
dimension of the energy-momentum tensor. For the bilocal operator T pu ' ,prT (x) the situation 
is more complicated. This is because of extra (quartic!) UV divergences appearing in the 
limit of x — >• 0. If x is kept away from then T pv ' ,pa (x) is finite and renormalization scheme 
independent. These divergences (which are local in x!) manifest themselves in the Fourier 
transform T p "' pa (q). They can and should be renormalized with proper counterterms: 

[T^{q)\ = T^{q) - ]T Z?{q)O h (18) 

i 

where 0{ are some operators of (mass) dimension < 4 and Zf t (q) are the corresponding 
(divergent) Z-factors. The latter must be local, that is have only polynomial dependence 
of the external momentum q. Within the MS-scheme Zf"(q) are just poles in e. It is of 
importance to note that the subtractive renormalization encoded in eq. (|18p is in general not 
constrained by the QCD charge renormalization. Thus, the unambiguous QCD predictions for 
the coefficient functions in OPE ([7]) could be made only modulo contact terms proportional 
to 5(x) in position space. 



3 OPE of the energy- momentum tensor correlator 

The leading coefficient Co is just the perturbative VEV of the correlator eq. ([T]) 

Cr pu {q) = (0|T^(g)|0)| pert (19) 

which we have computed up to order a 2 s (three loops) . In Figure [2] we show some sample 
Feynman diagrams contributing to this calculation. The energy-momentum tensor plays the 
role of an external current. In order to produce all possible Feynman diagrams we have used 



the program QGRAF 12]. As these diagrams are propagator-like the relevant integrals can 
be computed with the FORM package MINCER [18| after projecting them to scalar pieces. 
For the colour part of the diagrams the FORM package COLOR 19(] has been used. Because 
of the four independent external Lorentz indices there are many possible tensor structures 
for the correlator eq. (pQ) and hence the Wilson coefficients. These are composed of the large 
external momentum q and the metric tensor g. Using the symmetrical of eq. ([T]) we can 
narrow them down to five possible independent tensor structures: 

t pP°(q) = q p q"qPq°, 
tT Pa (q)=q 2 (Sfq v 9 p ° +q p q°9n, 

tT pa {q) = q 2 Wq p 9 va + + <z W 7 + gV<T) , (20) 



tT p \q) = (q 2 ) rr, 

tT' pa {q) = {q 2 )\9 pp 9 ua + g pu g vp )- 



These symmetries are fi < — > v,p < — > a and {pv) < — > (po")- 
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The conservation of the energy-momentum tensor leads to additional restrictions: 

q li T^ u ' f " T (q) = (local) contact terms (21) 

This condition has been checked in all calculations. Subtracting the physically irrelevant 
contact terms leads to only two independent tensor structures which have also been used e.g. 
in 



20]: 



t pu > pa {q) =if v r,P a 

t pU]pCT (q) =r fPrf a + r^rf" - -^—^rf v i m ^ 

with v P v {q)=q 2 gP u - q^q v 

where D is the dimension of the space time. The structure tjf ,pa (q) is traceless and orthogonal 
to tg U ' pa (q) . Hence the latter corresponds to the part coming from the traces of the energy- 
momentum tensors. The Wilson coefficients in eq. ([7]) are then of the general form 

r=1 ' B d to, (23) 

= C ;PCT («)(Q 2 )^^C| r) (Q 2 ) (+ contact terms). 

r=T,S 

where dim(Oo) = and dim(Oi) = 4 are the mass dimensions of the respective operators. 
The coefficients defined in the first line of eq. (|23p and there conversion to the ones defined 
in the second line are given in the appendix. In order to compute the coefficient C pu ' pa (q) 



we have used the method of projectors [2l|j which allows to express coefficient functions for 
any OPE of two operators in terms of massless propagator type diagrams only. For the case 
of C pu ' p<T (q) the following projector has been employed: 




C p J a {q) = P x \ III d 4 xdV d 4 y 2 e^+fcw+fcw {0\T[A«( yi ) A b v {y 2 ) f ^^(x)]|0) Q ~ irr 



(24) 
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Here 




J ki=0 



and the upper script Q~ irr means that only diagrams which become 1PI after formal gluing 
(depicted as a dotted line above) of two external lines carrying the (large) momentum q are 
included. Figure shows some sample diagrams at tree and one- loop level. 

a 




Figure 3: Diagrams for the calculation of the coefficient C\. 



4 Results 

a 2 

All results are given in the MS scheme with a a = a s = 4*- and the abbreviation l m = 

In (q ?) where ji is the MS renormalization scale. They can be retrieved from 

http: / /www-ttp. particle. uni- karlsruhe.de/Progdata/ttpl2/ttpl2-025/ 

The gauge group factors are defined in the usual way: C F and C A are the quadratic Casimir 
operators of the quark and the adjoint representation of the corresponding Lie algebra, d R 
is the dimension of the quark representation, n g is the number of gluons (dimension of the 
adjoint representation), T F is defined so that T F 5 ab = Tr (r a T b ^j is the trace of two group 

generators of the quark representation Q For QCD (colour gauge group SU(3)) we have 
C F = 4/3, C A = 3, T P = 1/2 and d R = 3. By n f we denote the number of active quark 
flavours. 

4.1 C 

Because of the contact terms both coefficients Cq and Cq could be unambiguously computed 
only up to constant (that is q-independent) contributions. To avoid the ambiguity we present 



4 For an SU(iV) gauge group these are d R = N,C A = 2T F N and C F = T F (N - i) . 
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below their Q 2 -derivatives: 



1 



dQ 



2^0 



16vr 2 

+ a a 



1 1 , 

n„ n f d R 

10 9 20 f 



1 



18 



A "-9 



1_ 
144' 



:n f T F n g 



67 



12960 
11 



2 3 
C> 9 + Y^ n f T FC F n g 



10663 ^ m 473 2rTl2 ^ 

n f C A T F n„ H n f T„n„ {^o) 

51840 / ^ 6480 



216^ 9 ^'' 4n£ ' 



1728 432 



11 3 1 

+ ^C3C> g + — ( 3 n f C A T F n g - — C 3 n 2 T|n 9 



^ dQ 2 w y 16vr 2 1 1296 A 3 162 5 9 81 ' J 



2 r7~~i2 

F^g 



144 7T 2 

where 



P0 n 9> 



(27) 



A - ~T2 3~ 

is the first coefficient of the perturbative expansion of the (3- function (|12|) . This result for 
Q 2 -^ji is in agreement with the one derived in [i^] for the case of gluodynamics (n/ = 0) 
at order g„ (two-loop level). The simple form of eq. (|27|) comes from the well-known trace 
anomaly [3, [z^] , which reads 

T% = ^ [Gp a G apa ] = -2/J(a.) [0 1 ]. (28) 
Indeed, from operator eq. ([25]) we expect that 

i J d 4 xe iqx (0|T[T^(x)T^(0)]|0) = 4 (3 2 {a s ) Q 4 TL GG (q 2 ) + contact terms, (29) 

where 

Q 4 U GG (q 2 )=i /"d 4 a;e i9a: (0|r[Oi(a;)Oi(0)]|0). (30) 
Now the one-loop result 




leads directly to eq. (|27|) The fact that this particular three-loop result can be derived from 
one- loop results is also the reason for the lack of (^-functions in it. Furthermore the structure 
of eq. (|28p explains nicely why the leading contribution for this scalar piece is of order a 2 . 
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In fact, the correlator (|30|) is known in two-, three- and four-loop approximations from works 
23], 24] and [l| respectively. The four- loop result reads (with all colour factors set to their 
QCD values and l m = 0) 



d 

d&< 



GGr„2\ 



H (q 



1 



16vr 2 



2 + a s 



73 7 
Y + 3^ 



+ at 



+ < 



+ n f 



37631 495 A 



15420961 44539 



864 

368203 11677 



-Cs 



7189 
~72~ 

3465 



-Cs 



108 



. 95 
2T C3 + 18 Cs 



115207 113 
648 + "llP 3 



7127 



127 2 
54 



2916 27 



(32) 



Finally, using eq. (|29|) and the well-known result for the four-loop QCD /3-function [25j, [26|] 
we could easily extend the rhs of (|27|) by three more orders in a s : 



Q 



+ a s 



d 



C { S \Q 2 ) 



121 22 



dQ 2 u ^ y 16tt 2 I 
( 11077 1025 



^ + ^ n/ "81^ 



+ 



27 
7 



72 36 
5787209 6655 



Tit -\ n*f 

162 / 243 y 



1728 



+ 



+ 



n 3 



16 
29071 



(3+rif 
5 



540049 4235 



.11664 162 
2351076745 5925007 



+ 



+ n f 
+ 



31104 288 
367411229 33359777 



648 72 
127 \ 



Ca 



+ nl 



4374 
46585 



16 



Cs 



15552 



240185 
7776 + 648 
381988321 3715127 , 

+ 7T7^ C3 



+ n 



/ 



139968 11664 
20279497 180083 . 95 



12485 
972 



-Cs 



139968 



17496 



Cs 



1458 



Cs 



1101389 427 
314928 + 2916 



7127 



556555 275 

1 

7776 108 s 



236196 2187 



Ca 



(33) 



4.2 Ci 



According to the definition of C^' P<T in eq. (|23|) there is a factor 7q2w in front of the dimen- 

(S) (T) 

sionless scalar pieces C\ and C\ which makes the whole coefficient immune to contact 
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terms except for those proportional to the tensor structures t^ u ' pa (q) and t pu,pa defined in 
eq. (|2U|) . The physical pieces C\ and C\ however are unambigous and the results read: 

°P = «■{ 

+ a2 .{§i C i-l n tT r Cr-^n,C A T r -^n 2 ,Tiy (34) 

+ a A-^ cl+ ^ n ' T ' c ' + ^ n ' c * T '-^' Ti \- (35) 

One thing to notice about C\ is that if we take the trace of both energy-momentum tensors 
the whole term i] p ^{q)f] u u {q) (Q-ip Cp^ in the Wilson coefficient becomes local and, therefore, 
indistinguishable from contact terms. We can however check eq. (|34j) independently by 
computing first the coefficient function C^ G,T ^ in an OPE ( tg U = q A , tjf = q 4 g^ u — 

i L 4 xe^T[T^(x)G 2 p M] — 

J q'— >— oo 



(36) 



(S) 



+ ... 



and then employing eq. (|28|) to get the next higher order in a s for C\ . The result 

,{TG,T) _ _ 16 (22 _ 8 \ 2 ^ _ _i6 fl ,„ \ /o\ _j_ 



Cr ' ; = -y + o s (^CU - -n f T F j + 0(af ) = -- (1 + P{a s )/2) + 0(af ) (37) 
allows to represent the rhs of eq. (|34|) in a form directly confirming eq. 



C f = fel C [ TG ' T) + 0{al) = ~ P(a.) (1 + P{a s )/2) + 0(a* ). (38) 

The factor /3(a s ) in this result is a direct consequence of the trace anomaly equation ([25]) . 
However, we do not know any rationale behind the peculiar structure after this factor. If it 
is not accidental, then one can hope that an explanation could be found within the so-called 



/3-expansion formalism suggested in [27]. 



(S) (T) 

It is important to note that the coefficient functions C\ and C\ are not Renormalization 
Group independent. We can construct the corresponding RG invariants by using the well- 
known faclo that the scale invariant version of the operator 0\ is 

(^' = ^a.)[(h\, Has) = Z ^ = a s \l + y j ^a\\. (39) 

From this and the scale invariance of T^ v ' ,( " T {q) defined in eq. ([T]) we find the RG invariant 
Wilson coefficients 

C\ S 1 CT = C{ s) /P(a s ) 

m m ( 4 °) 




5 This follows directly from the RG invariance of the energy-momentum tensor and the trace anomaly 
equation 
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which satisfy 

c[fl0^ = c[ S ' T \ 0l] . (41) 
Prom this definition we can immediately explain the absence of l m in eq. (|34p and eq. (|35p . 

(S) (S) x-rm 

Suppose we had l m in C\ and therefore in C\ then the general structure of eq. (|4"0j) up 
to three-loop order would be 



bl I ftg) + 0s(02 + &2 + C2 l^g) 

+ a 2 (a 3 + 63 l m + c 3 ^ 9 + d 3 1%) + O(a^) 



(42) 



with scale independent coefficients Oj,6j,Cj and di. The derivative with respect to fj 2 must 
vanish: 

^ 2 ^~2 C i S RGl = bi+ a s (b 2 + 2c 2 l m ) + a s /3(a s )(a 2 + 62 l m + c 2 Z^) 



d^ 2 

+ a 2 s (b 3 + 2c 3 i M + 3d 3 i^) + O(a^) = V^ 2 

61 = 

62 = 0, c 2 = 
=4> 63 = /3 a 2 , c 3 = 0, d 3 = 0. 



(43) 



In conclusion, not only have we explained the absence of logarithms in eq. (|34|) and eq. (|35|) 
but we also get the logarithmic part of the three-loop result for these coefficient functions 
for free. Terms with I 2 can only appear starting from four-loop level, terms with 1^ from 
five-loop level and so on. 

The quantities defined in eq. (|4U|) are given by 

Cfic, = §C A - ^n f T F - ^ (11C A - 4n f T F ) 2 (44) 

Ci% = + n M + 864(116^-4^) ( 2UC a + 876C 2 A n f T F 

+3537C A C F n f T F - 672C A n 2 T 2 - 1728C F n 2 T 2 + l&n^Tf) (45) 

The three-loop parts proportional to l m are 

r {S,2,l,\o & ) _ 2/ [ 1331C* 121 2 _ 11 2 2 _4_ a/rtfl ( 4fi \ 

^l,RGI ~ a s l M \ 3 888 ^ 324^ / A n f 1 F 81 L/AU f 1 F^2A3 n f 1 F(^ 

^(T,31,log) _ 2/ / 214C3+876C^n / T F +3537C jt C J rn / T F -672C A n^T|,-1728C F n^T|+16n3T| ) 
^l,RGI — "s'wj 10368 r-l 4 

For completeness we have also computed the contribution of the gluon condensate to the 
OPE of correlator ([BO]) : 

Q 4 n GG (g 2 ) = C7 GG Q 4 + Cf G (0|[Oi]|0) (48) 
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with the result: 

C GG 



l + a s 



"36 9 U ~ 12 m A %M n f lf 



+ a 
1151 



11509 2 13 3095 
97, 



-n 2 f T 2 F 
81 s 



-l m n i T F C F + -l m n f C A T F 



216 
121 



2y f f 



11 



33 



(49) 



144 



18 



-3( 3 n f T F C F + ^C3n/C^j,J 



+ 



The tree and one-loop contributions in (|49|) are in agreement with [8j and [28|, [29J] corre- 
spondingly. The two-loop part is new and has a feature that did not occur in lower orders, 
namely, a divergent contact term. Its appearance clearly demonstrates that non-logarithmic 
perturbative contributions to Cp G are not well defined in QCD, a fact seemingly ignored by 
the QCD sum rules practitioners (see, e.g. [j. liol]). 



An unambiguous QCD prediction can be made for the derivative: 

1151 



12 



1 T 

-n f T F 



216 



C 2 — rifTpCp — —n f C A T F + —n 2 T 2 



10 

27 5 



(50) 



~*~T2~^ 9 ^ -^-^g n jC A T F + -l M n 2 T 2 



5 Numerics 

In this section we will give our main results in the numerical form for two cases of interest, 
that is gluodynamics (nj = 0) and QCD with three light quarks only (nt = 3). As has already 
been mentioned, not all coefficient functions which we have discussed in the previous section 
are Renormalization Group independent. For a meaningful discussion we will construct the 
corresponding RG invariants by using the scale invariant version of the operator 0\ defined 
in eq. ([39]) . In addition we set l M = everywhere^ 

Q2 W 2 CP "80^2 i 1 ~ 1-66667 a s - 30.2162a 2 ) , (51) 

Q2 ~aW Ci ° S) ^Po "2^2°^ ( X + 22 - 8864as + 423 - 833a s + 8014.74af) , (53) 
Q 2 -^v? C o S) = -TTTTT^ f 1 + 18.3056a s + 247.48a 2 + 3386.41a^ , (54) 
6 This corresponds to the choice /i 2 = Q 2 for the renormalization scale. 



12 



r\1 d r<GG,RGI H 2/-, , m,™ s nGG.RGI _ 5/ \ nGG 

Q 17r2 C i == —a s (1 + 19.7576 a s ) , C x =f3(a s )C l , (55) 

C i GG ' KG/ — ; 7 a * C 1 + 15.3889 a s ) , (56) 

^iSg/=ST (i- 1 - 375 *-). C&i = ^ S V/3(a s ), (57) 

nf=\j y 

c S^^ 2(1 - L125as) ' (58) 

C'fie/^-^ 1 - - 2431825 ^)' C i%i = C l ] /ha s ), (59) 

— 3 (! " 1-3333 a,) . (60) 

6 Applications to high-temperature QCD 

Recently, the correlators II GG and T fJ,u,pa (q) have been studied in (Euclidean) hot Yang-Mills 



theory in 3l|, [32| respectively (see, also references therein for related earlier works). 



In this section we will employ our T = calculations in order to extend some of the results 
of these publications by adding fermionic contributions as well as higher order corrections. 
Note that for simplicity we will set all colour factors in all expressions below to their QCD 
values. The reader interested in expressions valid for generic colour group should be able to 
derive the corresponding results himself from our results. 

6.1 Trace anomaly correlator 

In 31] two- loop corrections to the quantitjQ 

G 9 (X) = (T[9(X)e(0)]) c , 6 = T% (61) 

where (. . . ) c stands for the connected part and the expectation value is taken at finite tem- 
perature^] T, have been computed. The capital case X for the space-time argument in (JBTj) is 
used in order to stress that we are dealing with a Euclidean correlator. In the following e and 
p = |q| are the energy and momentum densities with the well-known relation (9) c = e — 3p. 
In the limit of small r = \X\ the result of 31] reads H 



4a 2 „ , . 384 . , 8a s (9) c . . 64(e+p) , . ^/T 6 \ . . 



7 Note that Gg(0,X) has been directly measured in lattice simulations [33 |. 

8 We use the bold case for the temperature to make it distinct from T(. . . ) standing for the time ordered 
product of operators inside the round brackets. 

9 The expression below is the somewhat modified eq. (5.7) of [3l]]. 
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with 



79;i(r) = a 2 s +a 3 s f-^ + y + C(<# 

7 M (r) = 22a? + 0^), 
79;e+ P (r) = a^ + a^ + yl^ +0(a*), 



(63) 
(64) 
(65) 



and l„ x = log(^ 2 X 2 /4) + 2 lE 



According to [34| the coefficient functions jg-t and je;e( r ) do not depend on temperature T 
and, thus, should coincide with their T = counterparts. Hence, we can use our momentum 
space results described i n p revious sections to arrive at the following QCD predictions for 
junctions 10 !. 

%t(r) = a 2 + a* I --- + ^ l^ x + n f 



both coefficient funct 



1 11 

12 + T 



1-1/ 
18 3 mX 



49 495 



397 



"24 8 C3+ 16 



363 , 2 

16 >* X } 



35 5 43 11 2 

144 + 4 Cs " 12 mX " T > x 



+ n1 



216 + 36 ^ + 12 > x 



= / 255155 2915 A 3465 , 20891 , 
+ a\ ( — — —( 3 + — — C 5 + — / 



+ n f 



+ n ) 



+ n ) 



1728 
38741 9 



16 



192 



5445 

jiX ^ — C3 



1793 2 1331 3 



1728 16^ 36 ^ 5 



95 „ 16685 



576 



/ M x + 55C3 



4241 



— I 

96 mX 

361 125 491 5 289 2 11 3 

"216 + ~72^ + 1728 " X ~ 6 Cs ^ + 144 ^ x + 12 » x 

37 1 . 13_ _ J_ 2 _ J_, 
1458 27 ^ 3 + 324 Vx 108 ^ x " 



121 



54 



+ 0(a 



6 ] 

s ' - 



/ 3 



(66) 



7e ; e(0 = a 2 s 22 - - n f 



I 788 

+ (+— + 121Z M x+n/ 



304 44 



+ n 



4 

27 ' 9 



+ 0(aj). (67) 



Note that our vacuum calculations produce no information about the coefficient function 
70 ;e 4-p corresponding to the traceless part of the energy-momentum tensor. 

Numerically eqs. (fBUj) and (|57|) read (we set l^x = 0) 

(68) 



70.! (r) = aj -0.08333 < - 76.4189 o* + 82.4604a? + 0(a°), 

n f =0 

j e .t(r) = a 2 - 0.25 - 73.1821 + 142.705af + 0(a 6 s ), 

n f =3 



(69) 



%e(r) =22\a 2 + 11.9394 a? + 0(a% 



nt=0 



(70) 



The details of the corresponding Fourier transformation are spelled e.g. in [35 
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%e{r) — 18 [a 2 + 9.11111 a 3 s + 0(aj). (71) 



6.2 Shear stress correlator 



In 32|] the so-called shear stress correlator, defined as 

G V (X) = -16c 2 (T[T 12 (X)T 12 (0)]) c (72) 

with X = (Xq,X),X = (0,0, X3), has been computed up to two-loops in high-temperature 
Yang-Mills theory. Here is an arbitrary constant (introduced for some reason that is not 



quite clear to us in [32]]) which we put for simplicity equal to i/4. The calculation has been 
performed with the help of an ultraviolet expansion valid in the limit of small distances or 
large momenta; the result has been presented in the form of an OPE. As the corresponding 



Wilson coefficients should be T-independent the results of 32| can be checked and extended 
further with the help of our calculations 1^1 

We start from momentum space. In the zero temperature limit the function 

G V (Q 2 )= Jd 4 Xe^ x G v (X) 

is related to contribution to energy-momentum tensor correlator (pQ) proportional to the tensor 
structure t§ u,p(J {q). This fact could be easily checked by applying projector (2.5) of [32| to 
the correlator T )MU ' ,pa (q) expressed in terms of five independent tensor structures displayed in 
(1201). The result reads -8Q 4 (1 - 7/2 e + 7/2 e 2 - e 3 ) (cf (Q 2 ) + C|(Q 2 ) + ...). 

Thus, we will work with the representation 

T „ {q) t „ {g) { c s {q2) + C 5 (g2) (#) + \ + structures !_ 4 (73) 

q 2 ->-co y Q J 

We first concentrate on the coefficient function C|(Q 2 ) as the two- loop expression for C|(Q 2 ) 



presented in 32| is in agreement to the previously known expression obtained in [201 ] . The 



result of [32| for the second term in eq. ([7"3"|) reads: 

where Q12 is an unknown constant. Note that the second term of the above expression is 
obtained not from a calculation but with the use of Renormalization Group considerations 
similar to those leading to eq. (|43p . Such a derivation assumes that the coefficient function 
Cq is finite which is not obvious as the corresponding Feynman integrals have logarithmic 



n The Wilson coefficients in front of Lorentz non-invariant operators are for the moment not reachable with 
our projectors. It would be interesting however to extend these methods in order to reach e.g. the coefficient 
in front of (T 00 ) ~ e + p with a similar approach. 
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divergences stemming from the region of small x in eq. (UJ). Our direct calculation explicitly 
demonstrates the presence of such divergences: 



,(5)^2n 1 K / 41 7 \ 2 / 117 457 1 2 

lr(Q 2 ) = ; i + a 1 n , +a, n/H n, 

e w; 3/3 (a.) 1 V 24 288 V s 1 32 576 ; 576 7 



It is important to note that the contribution proportional to C| in (|73|) contains contact 
terms on/?/. This is in agreement with (|35|) due to an identity 

C (T) _ c (5) = CQntact termg) ( 76 ^ 

which, in turn, follows from restriction (|21|) (recall that Cg (Q ) = C\ /(— 2/3(a s )) as a 
consequence of ([25]) ). 

In Euclidean position space eq. (|73p can be presented as follows: 
fwiX) — ^ p 6™ + |c|(r) 1 + C|(r) <0|6>|0> + . . . j + structures 1-4 (77) 

Eq. (|76p . rewritten in terms of RG invariant quantities assumes the form: 

C i T RGi ~ 2 Po Cf> = contact terms. (78) 

By recalling that the contact terms do not contribute the function G^{x) for all x ^ we 
conclude that eqs. (|75p and ()47j> contain all information to construct the first non-zero term 
0(a 2 ) in the coefficient function Cq(x) with the result 

2 ftcf ( r ) = f H + 1Z n/ _ s 2 . i 3^ (79) 

MU 9W vr 2 r 4 ^192 16 ; 48 / 5184 f y v ; 

Finally, using the identity 

and (|26|) we arrive at the following result 



C P ~ C o 5) = contact terms, (80) 



^5, n 1 48 9 ( „ 7 \ 2/711 1188 

= ^8 y + 5 n / + G ^ " 16 +3 n / + G ? (-5 ^Cs-44/^ 



+ n f 



259 27 109 
120 ~ !) + ~12~ ^ X 



41 6 ^ 

90 + 5 C3 18^ 



(81) 



Numerical versions of eqs. ([79"]) and ([ST]) with Z M x = are presented below. 



,2 



107 



n/ =o tH r 4 192 



2ftC§(x) == _^^I = 0.557292}, (82) 



2/3 (5|(x) == -^(^=2.8125o), (83) 

n/=3 7T Z I 16 
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C\(x) = — -4-s \ 1 - 1-66667 a a - 14.9384a? V (84) 

n f =0 5 7T 4 r 8 y 3 I 

C\(x) = -3-5(1 -0.60.-10.69830?]. (85) 

n f =3 7T 4 r 8 \ I 



7 Discussion and Conclusions 

We have presented higher order corrections to coefficient functions Co and C\ of the OPE of 
two energy-momentum tensors in massless QCD as well as for the OPE of two scalar "gluon 
condensate" operators in massless QCD. Our results extend the previously known accuracy 
by one loop for the coefficient functions in front of the unit operator and by two loops for 
the CF of the gluon condensate operator 0\ = —\G tw G ilv . 

We have confirmed all previously available results and in some cases extended them from 
purely Yang-Mills theory to QCD. Contrary to previous assumptions, we have found that 
the coefficient functions Cf G as well as Cq 5 \q 2 ) are not completely finite with the standard 
QCD renormalization. 

We thank H. B. Meyer who has drawn our attention to the importance of the energy- 
momentum tensor correlator. Furthermore we would like to thank Y. Schroder, A. Vuorinen 
and M. Laine for useful comments. 

We are grateful to J. H. Kiihn for interesting discussions and support. 

In conclusion we want to mention that all our calculations have been performed on a SGI 
ALTIX 24-node IB-interconnected cluster of 8-cores Xeon computers using the thread-based 



361 ] version of FORM [37|] . The Feynman diagrams have been drawn with the Latex package 



Axodraw [3 

This work has been supported by the Deutsche Forschungsgemeinschaft in the Sonder- 
forschungsbereich/Transregio SFB/TR-9 "Computational Particle Physics". 

A Results for Cq^ and C[ 7 \ r = 1 . . . 5 and conversion to Cf )T) 

and C[ S,T) 

(r) (r) 

Here we give our intermediate results for the coefficients Cq and C\ (r = 1 . . . 5) appearing 
in the first line of eq. ([23 j) . i.e. the coefficients for the tensor structures eq. (|20|) before 
subtraction of contact terms: 

Ct ] =C? (P \ (86) 
Cf) =4 r C7f (r) . (87) 

The conservation of the energy-momentum tensor in classical field theory translates to 

T» v = local terms Q. (88) 

From this we get the relation 

1fiC^ v,pa {q) = (local) contact terms (i = 0, 1) (89) 
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which leads to the three restrictions serving as checks in our calculation: 



A (1) (Q 2 ) 
Df\Q 2 ) 



Ci,i(Q 2 ) + C ii2 (Q 2 ) + 2C it3 (Q 2 ) = (local) contact terms, 



(90) 



Ci^iQ 2 ) + Ci,4{Q 2 ) = (local) contact terms, 

Df ] {Q 2 ) := C ii3 (Q 2 ) +C ij5 {Q 2 ) = (local) contact terms. 

Hence the subtraction of contact terms enables us to write the Wilson coefficients in terms 
of only two independent tensor structures eq. (|22p which are related to the original five by 
the following equations: 



cf\Q 2 ) 
C\ T \Q 2 ) 



-CiMQ 2 ) - 



(D-l 



■C ii3 (Q 2 



(91) 



D^(Q 2 ) 
D$\Q 2 ) 



A.l Cj r) , r = 1...5 

These coefficients fulfill the relations eq. (|9Up even without local terms: 

D { q ] {Q 2 ) := C ,i(Q 2 ) + C ,2(Q 2 ) + 2C , 3 (Q 2 ) = 0, 

= CoMQ 2 ) + CoAQ 2 ) =0, (92) 
= C ,3(Q 2 ) + Co, 5 (Q 2 ) =0. 
Hence it is enough to give Co,4 and Co,5 here: 

16vr 2 ) C ,4 = +^ { _ i^44 a ^ n * + is^ '"' * 2 *""' ~ 38^8 a ^^ T ^ 9 } 

1/1 1 , 1 n 7 35 2 2 

H — < n„ n f a R H a,C A n„ a,n,l F n n a C,n„ 

e 1 15 9 30 ' 54 9 432 ' 9 11664 s 9 

1 2 tt n ^09 2 r< rr 77 2 2rp2 

H a n f l F C F n„ a n f C A l F n„ H a n,l F n„ 

192 s ' 9 93312 9 23328 

47 23 w 1 / 1 / w 

450 n9 225 %Q!fl 15^ 3o'«"' d * 
. J , 187 1987 1 7 

1 1620 9 12960 27 216 M9 7 9 



f 160831 2 61 1733639 140909 2 ^ 2 

l + 1399680 GAn9 + 480^°^ " 2799360 n ' C ^ n » + 699840^^^ 
941 n 2 1 ^ ^ 15943 593 2 ^2 

~*~9720 9 9 ~*~ 64 ni fUb ~ 77760 ^ n f A Fn s + iqUf ng 

1 r2 /-»2 1 109 ; 2 ^ rr-i _7_;2 2^2 

"648 A 9 + 5i84 f « n / GAij?Tl « " i296V«V* n « 

5 1 1 

-yfi$ C \ n B ~ lQ>n f T F C F n a + —(5n f C A T F n g 

563 A ^ 2 37 . m ^ 29 . „ m 19 2 ^2 

+ ^Cz^ q C 2 A n g + -^C3^ 9 re /CAT F n 9 - ^C3^ 9 ra 2 T|n s | , 

(93) 
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^167r 2N ) Cq 5 = H — k /h a 2 Cjn„ ^—a 2 n f C A T P n a H — a 2 n?THn l 

v y u > s e 2 \ 12 96 » ^ ff 10368 2592 s y F 9 J 

H — SH n„ H n f c? R a,C A n„ H a,n,T F n n a 2 C?n„ 

e I 10 9 20 ! 36 9 288 1 9 864 s A 9 

~ i^ ^ 7 ^™ 9 + 20^6 a ^ /CATFn9 ~ Slll ^^^ 9 } 
9 3 , 1 , 1 , 

+ + + j^ l ^ n 9 + ^ l ^ n f d R 

f 1 1367 1 7 

\ 540 9 8640 ' 9 18 M9 9 144 M9 7 9 

-^C3<? A n 9 - |^C3%r F n 9 | 

2 / 343429 2 307 983059 109129 2 ^ 2 

+ a < H C,n„ n f l F C F n„ H n f C A l F n n n f i „n„ 

s I 933120 A 9 1440 f 9 1866240 466560 / 9 

67 2 3 ™ n 10663 473 2t ,2 

a ^9 a 9 128 51840 6480 



12960 ™ A 9 i 2 8 w F ~ F " 9 ' 5 1840 /a-* f g 648Q < 
432^ C > B " 3456^ n/CATFns + 864' 



5 3 1 

gCsC 2 n 9 + -Csn/Tf-Cn,, - — ( 5 n f C A T F n g 

563 A ^ 2 3 ? > m ^ 29 A ^ m 19 2 ^2 

^Cs^a^ 9 - jgQC3W/r,C,n e + — C3n,CUT F n a + — £ 3n/ T F ) 



(94) 



A.2 C[ r) ,r = 1...5 

Here we give the five coefficents 



'4 7 

= a s { ~C A ~ J^ n f T F 



+ a l \^ n f T pC F + —n f C A T F - y^/^lj , 



(95) 



-C A + —UfT F 



n f 83 n 25 35 5 a al 

+ a 2 <^ C 2 + — TifTpCp + TiiC A T F + ?i 2 T 2 I , 

s \ 216 A 48 y 216 f 108 / F J ' 

+ a 2 ( — C 2 - —n f T F C F - —n f C A T F + — n 2 T 2 l , 
\ 432 96 f 432 ' 216 ; F J' 



(96) 



(97) 
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c 1A = + l 



— < — L A n,fl F 

e 1 36 9 1 

1 fl61 17 m I 

+ 3 + n2T6^-T08 n ^} ( 98 ) 

£ \ 72 12 ' 36 ' J 

+ a 2 i — C 2 n f T F C F n f C A T F n 2 T 2 l , 

3 1 16 A 144 7 108 y 108 1 F \ 



Ci,5 = - g 



7{~I5 - C7A + § a ' n ' T '} 



2 r 41 7 

3 +a ' {-108^-216^^ (") 

,2 



which fulfill the relations eq. (|90|) : 



£>f } =local / 0, (100) 
£>f } =local / 0. 

This is an important check as for the coefficent C^ u ' pa only counterterms of the form t^ u ' pa [Oi] 
and t^ v ' pa [Ox] are possible. Counterterms proportional to the other tensor structures would 
not be local. Hence = is necessary. 
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